Effective Gravitational Wave Stress-energy Tensor in Alternative Theories of Gravity 
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The inspiral of binary systems in vacuum is controlled by the stress-energy of gravitational ra- 
diation and any other propagating degrees of freedom. For gravitational waves, the dominant 
contribution is characterized by an effective stress-energy tensor at future null infinity. We employ 
perturbation theory and the short-wavelength approximation to compute this stress-energy tensor 
in a wide class of alternative theories. We find that this tensor is generally a modification of that 
first computed by Isaacson, where the corrections can dominate over the general relativistic term. In 
a wide class of theories, however, these corrections identically vanish at asymptotically flat, future, 
null infinity, reducing the stress-energy tensor to Isaacson's. We exemplify this phenomenon by 
first considering dynamical Chern-Simons modified gravity, which corrects the action via a scalar 
field and the contraction of the Riemann tensor and its dual. We then consider a wide class of 
theories with dynamical scalar fields coupled to higher-order curvature invariants, and show that 
the gravitational wave stress-energy tensor still reduces to Isaacson's. The calculations presented 
in this paper are crucial to perform systematic tests of such modified gravity theories through the 
orbital decay of binary pulsars or through gravitational wave observations. 

PACS numbers: 04.30.-w, 04.50.Kd, 04.25.Nx, 04.20.Fy 



I. INTRODUCTION 



Feynman has argued that no matter how beautiful or 
elegant a certain theory is, or how authoritative its pro- 
ponents, if it does not agree with experiments, then it 
must be wrong. For the past 40 years, this philosophy has 
been applied to gravitational theories with great success. 
Many modified gravity theories that were prominent in 
the 1970's, have now been essentially discarded, as they 
were found to disagree with Solar System experiments 
or binary pulsar observations Similarly, this decade 
is beginning to bring a wealth of astrophysical informa- 
tion that will be used to constrain new modified gravity 
theories. In fact, precision double binary pulsar observa- 
tions have already allowed us to constrain modified 
theories to exciting new levels 0, Q. Future gravita- 
tional wave (GW) observations on Earth, with the Ad- 
vanced Laser Interferometer Gravitational Observatory 
(aLIGO) 0-0, aVIRGO ^ and its collaborators, and in 
space, through the Laser Interferometer Space Antenna 
(LISA) [TlHS, ■^ill allow new precision tests of strong 
field gravity |15| . 

Such tests of alternative theories of gravity will be very 
sensitive to the motion of compact bodies in a regime 
of spacetime where gravitational fields and velocities are 
large, i.e. the so-called strong-field. Gralla [l^ has shown 
that motion in classical field theories that satisfy certain 
conditions (the existence of a Bianchi-like identity and 
field equations no higher than second-order) is "univer- 
sally" geodesic to leading-order in the binary system's 
mass-ratio, with possible deviations from geodesicity due 
to the bodies' internal structure. He also argues that one 
might be able to relax the second condition, as it does not 
seem necessary. In fact, motion in certain higher-order 
theories, such as Chern-Simons modified gravity [l^l, is 



already known to be purely geodesic to leading-order in 
the mass-ratio, without infiuence of internal structure 
due to additional symmetries in the theory. 

Tests of modified gravity theories in the strong-field, 
however, not only require a prescription for the conser- 
vative sector of motion, but also of the dissipative sector, 
that which describes how the objects inspiral. Geodesic 
motion must thus be naturally corrected by a radiation- 
reaction force that drives non-geodesic motion toward an 
ultimate plunge and merger jlSj . Similarly, one can think 
of such motion as geodesic, but with varying orbital el- 
ements (T9l - [23j (energy, angular momentum and Carter 
constant). The rate of change of such orbital elements 
is governed by the rate at which all degrees of freedom 
(gravitational and non-gravitational) radiate. 

In the gravitational sector and to leading order in the 
metric perturbation, such a rate of change is controlled 
by an effective stress-energy tensor for GWs, first com- 
puted by Isaacson in General Relativity (GR) [13, HH. In 
his approach, Isaacson expanded the Einstein equations 
to second order in the metric perturbation about an arbi- 
trary background. The first-order equations describe the 
evolution of gravitational radiation. The second-order 
equation serves as a source to the zeroth-order field equa- 
tions, just like a stress-energy tensor, and it depends on 
the square of the first-order perturbation. This tensor 
can then be averaged over several gravitational wave- 
lengths, assuming the background length scale is much 
longer than the GW wavelength (the short-wavelength 
approximation). In this approximation, Isaacson found 
that the effective GW energy-momentum tensor is pro- 
portional to the square of first partial derivatives of the 
metric perturbation, i.e. proportional to the square of 
the gravitational frequency. Components of this stress- 
energy then provide the rate of change of orbital ele- 
ments, leading to the well-known quadrupole formula. 
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Alternative theories of gravity generically lead to a 
modified effective GW stress-energy tensor. It is some- 
times assumed that this stress-energy tensor will take the 
same form as in GR [l^H^I, but this need not be the case. 
In GR, the scaling of this tensor with the GW frequency 
squared can be traced to the Einstein-Hilbert action's 
dependence on second-derivatives of the metric pertur- 
bation through the Ricci scalar. If the action is modified 
through the introduction of higher-powers of the curva- 
ture tensor, then the stress-energy tensor will be pro- 
portional to higher powers of the frequency. Therefore, 
the consistent calculation of the modified Isaacson ten- 
sor needs to be carried out until terms similar to the GR 
contribution (proportional to frequency squared) are ob- 
tained. This in turn implies that calculations of effective 
energy-momentum tensors in modified gravity theories 
to leading- order in the GW frequency can sometimes be 
insufficient for determining the rate of change of orbital 
elements. 

In this paper, we present a formalism to compute the 
energy-momentum tensor consistently in generic classi- 
cal field theories. We employ a scheme where the ac- 
tion itself is first expanded in the metric perturbation 
to second order, and the background metric and metric 
perturbation are treated as independent fields. Varying 
with respect to the background metric leads to an effec- 
tive GW stress-energy tensor that can then be averaged 
over several wavelengths. This produces results equiva- 
lent to Isaacson's calculation. 

We exemplify this formulation by first considering CS 
gravity This theory modifies the Einstein-Hilbert 

action through the addition of the product of a scalar field 
with the contraction of the Riemann tensor and its dual. 
This scalar field is also given dynamics through a kinetic 
term in the action. The leading-order contribution to the 
CS-modified GW stress-energy tensor should appear at 
order frequency to the fourth-power, but Sopuerta and 
Yunes |28| have shown that this contribution vanishes at 
future null infinity. 

We here continue this calculation through order fre- 
quency cubed and frequency squared and find that such 
CS modifications still vanish at future null infinity. This 
is because the background scalar field must decay at a 
certain rate for it to have a finite amount of energy in 
an asymptotically-flat spacetime. If one insists on ig- 
noring such a requirement, such as in the case of the 
non-dynamical theory, then frequency-cubed CS modifi- 
cations to the energy-momentum do not vanish. 

We explicitly calculate such modifications for a canon- 
ical embedding, where the scalar field is a linear function 
of time in inertial coordinates. This is similar to pre- 
vious work [2^ that calculated another effective stress- 
energy tensor for the non-dynamical version of Chern- 
Simons. In this case, the dominant modification to the 
radiation-reaction force is in the rate of change of radi- 
ated momentum, which leads to so-called recoil veloci- 
ties after binary coalescence. In GR, such recoil is pro- 
portional to the product of the (mass) quadrupole and 



octopole when multipolarly decomposing the radiation 
field. In non-dynamical GS gravity with a canonical em- 
bedding, the recoil is proportional to the square of the 
mass quadrupole, which dominates over the GR term. 

We then construct a wide class of alternative theories 
that differ from GR through higher order curvature terms 
in the action coupled to a scalar field. We compute the 
GW stress-energy-momentum tensor in such theories and 
find that corrections to the Isaacson tensor vanish at fu- 
ture null infinity provided the following conditions are 
satisfied: (i) the curvature invariants in the modification 
arc quadratic or higher order; (ii) the non-minimally cou- 
pled scalar field is dynamical; (iii) the modification may 
be modeled as a weak deformation away from GR; (iv) 
the spacetime is asymptotically flat at future null infin- 
ity. These results prove that the effective stress-energy 
tensor assumed in |26l . [27| is indeed corrcclQ. 

Even if the effective GW energy-momentum tensor is 
identical to that in GR, in terms of contractions of first 
derivatives of the metric perturbation, this does not im- 
ply that GWs will not be modified. First, background 
solutions could be modified. For example, in dynamical 
GS gravity, the Kerr metric is not a solution to the mod- 
ified field equations for a rotating black ho le ( BH) [30| . 
but it is instead modified in the shift sector [3l|. Second, 
the solution to the GW evolution equation could also 
be modified. For example, in non-dynamical GS grav- 
ity, GWs become amplitude birefringent as they propa- 
gate [s^-dll- Third, additional degrees of freedom may 
also be present and radiate, thus changing the orbital 
evolution. All of these facts imply that even if the Isaac- 
son tensor correctly describes the effective GW energy- 
momentum tensor, GWs themselves can and generically 
will be modified in such alternative theories. 

In the remainder of this paper we use the following 
conventions. Background quantities are always denoted 
with an overhead bar, while perturbed quantities of first- 
order with an overhead tilde. We employ decompositions 
of the type g^^ = g^^ -f- e h^^ + 0{e^), where g^^ is 

the full metric, g^^ is the background metric and hf^, 
is a small perturbation (e <C 1 is a book-keeping pa- 
rameter). Covariant differentiation with respect to the 
background metric is denoted via V^-B^, while covariant 
differentiation with respect to the full metric is denoted 
via V^i?^. Symmetrization and antisymmetrization are 
denoted with parentheses and square brackets around the 
indices respectively, such as ^(^j,) = [^^1/ + ^,y^]/2 and 
Aj^^j = [A^^ — Ajj^]/2. We use the metric signature 
(— , -|-, -|-, -|-) and geometric units, such that G = c = 1. 

This paper is organized as follows: Section lU describes 
the perturbed Lagrangian approach used in this paper 
to compute the effective GW stress-energy tensor. Sec- 



Tho authors of |26| . I27I presented an energy loss formula which 
was not evaluated at . In the limit of r — > 00, their energy 
loss formula reduces to the Isaacson formula. 
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tion mil applies this framework to GR. Section IIVI dis- 
cusses dynamical CS gravity. Section|V]computes the full 
effective stress-energy tensor in this theory. Section IVII 
generalizes the calculation to a wider class of alternative 
theories. Section IVlIl concludes and points to future re- 
search. 

11. PERTURBED LAGRANGIAN APPROACH 

Isaacson [23, HH introduced what is now the stan- 
dard technique to obtain an effective stress-energy tensor 
for gravitational radiation, via second-order perturbation 
theory on the equations of motion. This technique re- 
quires an averaging procedure to construct an effective 
stress-energy tensor. This is because of the inability to 
localize the energy of the gravitational field to less than 
several wavelengths of the radiation, and because of the 
ambiguities of the metric perturbation on distances of or- 
der the wavelength due to g auge freedonjf]- Isaacson em- 
ployed the Brill-Hartle [33] averaging scheme, although 
one can arrive at an identical quantity by using different 
schemes (ssl - lioj . e.g. Whitham or macroscopic gravity. 

An alternative approach to derive field equations and 
an effective stress-energy tensor for GWs is to work at the 
level of the action. One possibility is to use the Palatini 
framework [4l|, \^ , where the connection is promoted to 
an independent field that is varied in the action, together 
with the metric tensor. Such a framework, however is 
problematic in alternative theories of gravity, as it need 
not lead to the same field equations as variation of the 
action with respect to the metric tensor only. 

A similar but more appropriate approach is that of 
second-variation [4^ . In this approach, the action is first 
expanded to second order in the metric perturbation, as- 
suming the connection is the Christoffel one T'^^. Then, 
the action is promoted to an effective one, by treating the 
background metric tensor and the metric perturbation as 
independent fields. Variation of this effective action with 
respect to the metric perturbation and the background 
metric yields the equations of motion. The former vari- 
ation leads to the first-order field equations, when the 
background field equations are imposed. The latter vari- 
ation leads to the background field equations to zeroth 
order in the metric perturbation and to an effective GW 
stress-energy tensor to second order. 

Let us begin by expanding all quantities in a power 
series about a background solution 

= (^ + e(^ + e^<^ + 0(e^) , (1) 

where e ^ 1 is an order counting parameter and ip rep- 
resents all tensor fields of the theory with indices sup- 
pressed: If is the background field, is the first-order 



^ Gauge freedom in perturbation theory stands for the freedom 
to identify points between the physical and "background" man- 
ifolds [Hlsl. 



perturbation to (/S, and (f is the second-order perturba- 
tion. The action can then be expanded, as 

s[if] ^ m + ^(1) [^, ^] + 5(2) ^] + o(,3) ^ 

where S"*-^' is linear in ip and 5*'^^ is quadratic in (p but 
linear in (p. We now define the effective action as Eq. ([2]) 
but promoting (p and (p to independent fields. 

One might wonder why the field ip is not also treated 
as independent. First, variation of the action with re- 
spect to (p would lead to second-order equations of mo- 
tion, which we are not interested in here. Second, the 
variation of the action with respect to (p cannot intro- 
duce terms that depend on ip, because the product of 
ip and ip> never appears in Eq. ([2]), as this would be of 
0{e^). Third, the variation of the action with respect 
to ip can only introduce terms linear in (p, which vanish 
upon averaging, as we describe in Sec. Ill Al This is be- 
cause averages of any odd number of short-wavelength 
quantities generically vanish. Therefore, we can safely 
neglect all terms that depend on (p in the effective ac- 
tion, which renders Eq. ([2]) a functional of only (p and 
>p>. 

As in the standard approach, the second-order vari- 
ation method still requires that one performs a short- 
wavelength average of the effective stress-energy tensor. 
Upon averaging, the variation of the first-order piece of 
the action S^^-* with respect to (p> vanishes because it gen- 
erates terms linear in ip. Since S^^^ does not contribute 
to the effective stress-energy tensor, we can safely drop 
it from the effective action for now. 

The effective action reduces to 

5-[^,^]=5-W[^] + 5-(2)[^,^]. (3) 

Naturally, the variation of S*""'"^ with respect to the back- 
ground metric g'^" yields the background equations of 
motion. The effective stress-energy tensor comes from 
averaging the variation of S''=''(^) with respect to g^'^ , 

55°"^'^ J d'xV^ Sr^t^, , (4a) 
T-^-262((V)) , (4b) 

where the factor of 2 is conventional for agreement with 
the canonical stress-energy tensor, and (( )) is the av- 
eraging operator, which we discuss below. One of the 
immediate benefits of working from an action principle 
comes from the diffeomorphism invariancc of the action. 
The diffeomorphism invariance immediately implies that 
the variation of the total action with respect to the met- 
ric is divergence free [13] . When the matter stress-energy 
tensor is itself divergence free, then the gravity sector ~ 
the sum of the stress-energy of non-minimally coupled 
degrees of freedom and the effective stress-energy tensor 
of gravitational waves - will also be divergence free. 
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A. Short-wavelength averaging 

The goal of the averaging scheme is to distinguish ra- 
diative quantities, those which are rapidly-varying func- 
tions of spacetime, from Coulomb-like quantities, those 
which are slowly- varying functions. This is accomplished 
by defining the operator (( )) as a linear integral operator. 
This operator may either be an average over the phase of 
the rapidly varying quantities or over spacetime. If the 
integral is over spacetime, there is an averaging kernel 
with characteristic length scale Lavo that separates the 
foreground short-wavelength Aqw from the background 
length scale ibg, that is, Agw ^ -^ave < L^g- 

The details of the avera ging sc heme are not as im- 
portant as their properties |38l44g |. since one arrives at 
equivalent results using different schemes. The most use- 
ful properties of (( )) are: 

1. The average of a product of an odd number of 
short-wavelength quantities vanishes. 

2. The average of a derivative of a tensor vanishes, e.g. 
for some tensor expression T^^^, ((V^T'^q^)) = 0. 

3. As a corollary to the above, integration by parts can 
be performed, e.g. for tensor expressions R^a,Si3, 

Let us briefly mention where some of these proper- 
ties come from and some caveats. When considering 
monochromatic functions, an odd oscillatory integral has 
an average value about zero, while an even oscillatory in- 
tegral has a non-zero average. This is enough to find that 
that averages of expressions linear in a short- wavelength 
quantity will vanish. Expressions at third (and higher 
odd) order would vanish for monochromatic radiation, 
but not in general. However, these are at sufficiently 
high order that we neglect them. 

The vanishing of averages of derivatives is a subtle 
point. In the spacetime average approach, this is found 
by integrating by parts, leaving a term with a derivative 
on the averaging kernel. This term is smaller than non- 
vanishing averages by a factor of ©(Acw/iavc) and de- 
pends on the averaging kernel. From physical grounds, 
the choice of averaging kernel should not affect any phys- 
ical quantities, so the average should in fact vanish iden- 
tically. From the action standpoint, the average of a 
derivative can be seen to arise from an action term which 
is a total divergence. Since total divergences in the ac- 
tion do not affect the equations of motion, the average of 
a derivative vanishes. 

A similar argument holds for integration by parts. In 
the Brill-Hartle average scheme, integration by parts in- 
curs an error of order ©(Acw/iave) from a derivative 
of the averaging kernel. From the action standpoint, 
though, integration by parts at the level of the action 
incurs no error, since there is no averaging kernel in the 
action. 

The fact that the variation of S^^^ with respect to (p 
does not contribute to the effective stress-energy tensor is 



a direct consequence of property (1) above. The vanish- 
ing of all terms linear in if upon averaging is also a con- 
sequence of (1). As one can see, these properties greatly 
simplify all further calculations. 



B. Varying Christoffel and curvature tensors 

Let us now consider what types of terms arise from 
the variation of the effective action with respect to g'^'^ . 
In order to perform this variation properly, any implicit 
dependence of the action on g^^" must be explicitly re- 
vealed; for example, terms which contain the trace h must 
be rewritten as g^i,h^'^ . Indices should appear in their 
"natural" positions (see Sec. Ill D|) . for indices raised and 
lowered with the metric have implicit dependence on it. 
Furthermore, since we are approaching gravity from the 
metric formulation, rather than the Palatini formulation, 
there will be contributions from the variation of and 
curvature tensors. 

Consider a term in the effective action such as 



(5) 



where g is the determinant of the metric and 
T'^'" s---, S"'"' p... are some tensor expressions, with all in- 
dices contracted in some fashion, as the action must 
be a scalar. When such a term is varied with g'"^ 
g'^^ +5g'^'^, besides the obvious contributions from 
and explicit dependence of T'^"'s..., S"-'" p... on g''", there 
are also contributions from varying the Christoffel con- 
nection in Vp. The general expression is 

<5(V^5"--""/3,..;3„) - V^(<5^"^-">,.../3„)+ 

n rn 



where • . . A • • • in the i"^ term of a sum means replacing 
a; or Pi in the index list with A; and where 



ST" = 



=Q/3 



(7) 

Curvature tensors also depend on derivatives of the 
connection, so one naturally expects terms of the form 
VpVCTfSg^" from the variation of curvature quantities, 
i.e. the Riemann tensor R^^ajSi Ricci tensor i?^^, or Ricci 
scalar R. For example, one can show that 



(8) 



where the contribution from F A F cancels [42]. Upon 
integration by parts any scalar in the action that con- 
tains curvature tensors, one can convert a term contain- 
ing VpVo-fJg'^" into 



-g P^.^aSg'"' . 



(9) 
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In fact, many terms in the variation of the action can be 
written in the form of Eq. ^ . 

The contribution of Eq. ([9]) to the effective stress- 
energy tensor is found by integrating by parts and then 
averaging, according to Eq. (|4]). Upon averaging, how- 
ever, one finds that such terms do not contribute to the 
effective stress-energy tensor because 

t;«-- = 2((v.p%,)) , (10) 

vanishes according to property (2) in Sec. Ill Al 

The above arguments and results imply that the vari- 
ations of curvature tensors and connection coefficients 
with respect to g^,^ do not contribute to the effective 
stress-energy tensor. Only metric tensors which are rais- 
ing, lowering, and contracting indices in the action con- 
tribute to this tensor. We can thus concentrate on these, 
when computing T°". 

C. Contributions at asymptotic infinity 

When calculating the radiation-reaction force to lead- 
ing order in the metric perturbation, it is crucial to 
account for all the energy-momentum loss in the sys- 
tem. The first contribution is straightforward: energy- 
momentum is radiated outward, toward future, null in- 
finity, J^+. Since the stress-energy tensor is covariantly 
conserved, the energy-momentum radiated to J^~^ can be 
calculated by performing a surface integral over a two- 
sphere at future, null infinit}|f|. 

However, not all energy-momentum loss escapes to 
infinity, as energy can also be lost due to the pres- 
ence of trapped surfaces in the interior of the spacetime. 
Trapped surfaces can effectively absorb GW energy- 
momentum, which must also be accounted for, e.g. in 
the calculation of EMRI orbits around supermassive 
BHs [13, S^. Calculations of such energy- momentum 
loss at the BH horizon are dramatically more compli- 
cated than those at and we do not consider them 
here. 

What is the relative importance of energy-momentum 
lost to ^+ and that lost into trapped surfaces? To an- 
swer this question, we can concentrate on the magnitude 
of the leading-order energy flux, as the argument trivially 
extends to momentum. The post-Newtonian (PN) ap- 
proximation (46j , which assumes weak-gravitational fields 
and slow velocities, predicts that the energy fiux carried 
out to is proportional to u^" to leading-order in v, 
where v is the orbital velocity of a binary system in a 
quasi-circular orbit (see e.g. [4a]). On the other hand, a 
combination of the PN approximation and BH perturba- 
tion theory predicts that, to leading-order in v, the en- 
ergy flux carried into trapped surfaces is proportional to 



We will not consider spacetimes which are not asymptotically 
flat, e.g. de Sitter space; the calculations are more involved in 
such spacetimes. 



v^^ for spinning BHs and v^^ for non-spinning BHs (47| . 
BH GW flux absorption is then clearly smaller than the 
GW flux carried out to if u < 1, which is true for 
EMRIs for which the PN approximation holds. 

Intuitively, this hierarchy in the magnitude of energy- 
momentum flux lost by BH binaries can be understood 
by considering the BH as a geometric absorber in the 
radiation field. Radiation which is longer in wavelength 
than the size of the BH is very weakly absorbed. Only 
at the end of an inspiral will the orbital frequency be 
high enough that GWs will be significantly absorbed by 
the horizon. Notice that this argument is independent 
of the particular theory considered, only relying on the 
existence of trapped surfaces. This result docs not imply 
that BH absorption should be neglected in EMRI mod- 
eling, but just that it is a smaller effect than the flux 
carried out to infinity [3, . 

In the remainder of the paper, wc will only address 
energy-momentum radiated to ^+ and relegate any anal- 
ysis of radiation lost into trapped surface to future work. 
The only terms which can contribute to an energy- 
momentum flux integral on a 2-sphere at are those 
which decay as r~'^, since the area element of the sphere 
grows as r^. No terms may decay more slowly than r~^, 
as the flux must be finite, i.e. the effective stress-energy 
cannot scale as r"-'^, as a constant or with positive pow- 
ers of radius. Similarly, any terms decaying faster than 
do not contribute, as they would vanish at J^~^. Of 
course, to determine which terms contribute and which 
do not, one must know the leading asymptotic forms of 
all quantities in the effective stress-energy tensor. 

In GR, as we shall see in ISec. III[ the only fields ap- 
pearing in the effective stress-energy tensor are the back- 
ground metric g^j^ and derivatives of the metric perturba- 
tion hf^i,. As one approaches , g^^ ~ ry^i/ in Cartesian 
coordinates, while ~ ~ |Vp/).^,y|. Curvature 

tensors scale as \Rfj.^S(7\ ^ since they quantify tidal 
forces. For a theory that is a deformation away from 
GR, and far away from regions of strong curvature, these 
asymptotic forms cannot change. 

Consider now terms in the general effective action at 
order O(e^) that contain background curvature tensors. 
Due to their ordering, they would contribute to the ef- 
fective stress-energy. One such term is 



= / d^xV^ Vp^i Va(^f R'p.^ , (11) 



where <^i,<p2 are the first-order perturbations to two 
fields in the theory (e.g. the metric perturbation h^i, 
and the CS scalar perturbation 'd that we introduce in 
ISec. V|) . Since there is no contribution to the effective 
stress-energy tensor from the variation of curvature quan- 
tities f see ISec. II B[) . the only contributions to the effec- 
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tive stress-energy comes from 



t; 



1 



Vp(^? V„<^f r'p.^ ) ) , (12) 



which has the same functional form as the integral. Note 
that the curvature tensor always remains when varying 
with respect to g'"^. 

Combining this result with the asymptotic arguments 
above, such terms can be ignored as one approaches J^^ . 
Each of the first-order fields possess a radiative part that 
scales as r~^. The square of the first-order fields would 
then satisfy the scaling requirement for the flux inte- 
gral. The curvature tensor, however, scales as r~'^, which 
implies that the term in Eq. vanishes at 

We then conclude that terms in the action that contain 
background curvature quantities at 0{e^) may be ignored 
in calculating the effective stress-energy tensor at 
As an immediate corollary to this simplification, we may 
also freely commute background covariant derivatives if 
we are interested in the stress-energy tensor at infinity 
only, since the commutator is proportional to background 
curvature tensors. 



D. Imposing gauge in the effective action 

We will choose as our dynamical field not hf^^, but 
rather h , where the underline stands for the trace- 
reverse operation, and we take the "natural" position 
of the indices to be contravariant. The resulting stress- 
energy tensor is equal to the one calculated using h^i, 
after evaluating both of them on-shell, i.e. imposing the 
equations of motion. 

We also impose a gauge condition to simplify future 
expressions: the Lorenz gauge condition. 



V7 If^" 



0. 



(13) 



Typically one may not impose a gauge condition at the 
level of the action. However, in our case, the gauge condi- 
tion in Eq. (jl3p has the important property of having all 
of the indices in their natural positions: the contraction 
of the indices does not involve the metric. 

Consider a term in the effective action that contains 
this divergence. 



dS 



-g Tp v„r^ . 



(14) 



with Tp some tensor expression at first order in e. The 
a index that is contracted above does not require the 

metric for such contraction. Therefore, Vq/i always 
remains upon variation. 



rricff.c 



-2e' 



5T0 \ ~ »(3 



(15) 



If we delayed imposing the Lorenz gauge condition un- 
til after the calculation of the effective stress-energy ten- 
sor, we would find the same effective tensor as if we had 
imposed the gauge condition at the level of the action. 
Having said that, one should not impose the gauge condi- 
tion when varying with respect to as clearly Vah ^ 
must also be varied. 



III. EFFECTIVE STRESS-ENERGY IN GR 

Let us now demonstrate the principles described in the 
previous section by deriving the standard Isaacson stress- 
energy tensor in GR. Consider the Einstein-Hilbert ac- 
tion, 



'5'gr = d**2;y^ R, 



(16) 



where k = (IGttG) ^. Now perturb to second order to 
form the effective action, 



Qctt _ c"="'(0) I c"=ff(2) 



Q<=tt(0) 
•Jgr 



<-,ctt(2) 



(17a) 

(17b) 
(17c) 



where 



and 



4i?, 



al3 (^h ^jl 

( 



h n 



(17d) 



-F^V^V^F^ - i(V^^)(V^^)- 



(VpFa)(vX") + ^(v.F(v^F') -F^VpV^F^ 



, (17e) 



and {R^ij^R) refer to the background Ricci tensor and 

scalar respectively. The integrands have been written 

~ pi/ 

in terms of the trace-reversed metric perturbation, h . 
From lSec. II C[ £gr^ does not contribute at -J^^ because 
it depends explicitly on curvature quantities, so we ignore 
it. The variation and averaging of £gr^ produces the 
Isaacson stress-energy tensor. 

By integrating by parts, all terms in Eq. (|17e|) can 

be written as (VaF ){^ ^) (with indices contracted 

to form a scalar) rather than VaV/jft- (again, with 
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indices contracted). £gr is thus rewritten in the more 
compact form 



4(VmA„,)(vT'') + ^(v4)(v^A) 



(18) 



which is the expression that appears in MacCahum and 
Taub [4^. With this simphfied expression at hand, we 
can promote h to an independent dynamical field in 
Eq. ([T7)) and vary it with respect to both g^^ and ji 
to obtain the effective stress-energy tensor and the first- 
order equations of motion respectively. 

Let us first derive the first-order equations of motion. 
Varying Eq. (flT)) with respect to ji , we find 



Let us now evaluate the MacCallum-Taub ten- 
sor on shell, by imposing the equations of motion 
[Eq. (fT9|) ]. When short- wavelength averaging, deriva- 
tives that are contracted together can be converted into 
the d'Alembertian via integration by parts; such terms 
vanish at J^"''. What results is the usual Isaacson stress- 
energy tensor. 



GR/J,Z^ ^ 



(20c) 



Notice that this expression is only valid at and in 
TT gauge. As mentioned earlier, this tensor cannot be 
used to model energy-momentum loss through trapped 
surfaces, since then curvature quantities cannot be ig- 
nored. 



□Zi^. - 2V V(^Zi,)„ 



1 



(19a) 



IV. CHERN-SIMONS GRAVITY 



whose trace is 



(19b) 



We can now impose the Lorenz gauge on Eq. (|T9b)) . which 
then leads to Q/i = 0. If /i = is further imposed on an 
initial hypersurface while maintaining Lorenz gauge, then 
the evolution equation preserves the trace- free gauge [4l| . 
The combination of these two gauge choices (Lorenz 
gauge plus trace-free) is the transverse-tracefree gauge, 
or TT gauge. 

After commuting derivatives in Eq. (|19ap and imposing 
TT gauge, the tensor equation of motion reads 



(19c) 



where R^ai^p is the background Riemann tensor. At 
this equation reduces to Q/i^j, = 0, which leads to the 
standard dispersion relation for GWs, traveling at the 
speed of light. 

Let us now calculate the effective stress-energy tensor. 
Note that the first term in may be integrated by 
parts and covariant derivatives commuted to form the 
Lorenz gauge condition, so the first term may be ignored. 
Varying the action with respect to g^"^, we find 



where 



1 



1 — ~a/3 — - 



— n'^m-Tiiu 1 (20a) 



1-. 



+ ^^^^^ 



(20b) 



which we refer to as the MacCallum-Taub tensor. Terms 
that depend on the trace /i ^ in this tensor can be elim- 
inated in TT gauge. 



CS gravity is a modified theory introduced first by 
Jackiw and Pi [s^ (for a recent review see [l3l)- The 
dynamical version of this theory modifies the Einstein- 
Hilbert action through the addition of the following 
terms: 



S = S 



where 



Q — 



EH 



CS 



(21) 



:^ I d^xV^ ^ *RR, 



(V^^9) (V,i9) , 



d*x £n 



(22) 



The quantity k — (IGttG)^^ is the gravitational constant, 
while a and /? are coupling constants that control the 
strength of the CS coupling to the gravitational sector 
and its kinetic energy respectively. In the non-dynamical 
version of the theory, /3 = and there are no dynamics for 
the scalar field, which is promoted to a prior-geometric 
quantity. 

The quantity d is the CS field, which couples to the 
gravitational sector via the parity-violating Pontryagin 
density, *i? i?, which is given by 

ie"^^''i?„^,,i?''^^,, (23) 

where the asterisk denotes the dual tensor, which we con- 
struct using the antisymmetric Levi-Civita tensor s""^^^ . 
This scalar is a topological invariant, as it can be written 
as the divergence of a current 



UK .— -K^p^g H ^ ' — —e 



*RR = 4V, 



1 ^ 

7cr ^2 ^^V 7^ 



(24) 
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Equation ([21]) contains several terms that we describe 
below: the first one is the Einstein-Hilbert action; the 
second one is the CS couphng to the gravitational sector; 
the third one is the CS kinetic term; and the fourth one 
stands for additional matter degrees of freedom. The 
CS kinetic term is precisely the one that distinguishes 
the non-dynamical and the dynamical theory. In the for- 
mer, the scalar field is a priori prescribed, while in the 
dynamical theory, the scalar field satisfies an evolution 
equation. 

The field equations of this theory are obtained by vary- 
ing the action with respect to all degrees of freedom: 



C + —C = — ( T""" 



a 



(25a) 
(25b) 

where T^J"* is the matter stress-energy tensor and T^V is 
the CS scalar stress-energy: 



T. 



(V^i9)(V,i?) 



,,(V'^i5)(V.t?) 



(26) 



The C-tensor C^'^ is given by 

= {V J) e''^''^°'V ^R^\ + (V^V^i?) *i?«("'9)'^ . (27) 

Many solutions to these field equations have been 
found. In their pioneering work, Jackiw and Pi 
showed that the Scwharzschild metric is also a solu- 
tion in CS gravity [32j . Later on, a detailed analy- 
sis showed that all spherically symmetric spaces, such 
as the Friedman-Robertson- Walker metric, are also so- 
lutions [3^. Axially-symmetric spaces, however, are 
not necessarily solutions, because the Pontryagin den- 
sity does not vanish in this case, sourcing a non-trivial 
scalar field. Specifically, this implies the Kerr metric is 
not a solution. 

A slowly-rotating solution, however, does exist in dy- 
namical CS gravity. Yunes and Pretorius [3l| found that 
when the field equations are expanded in the Kerr pa- 
rameter a/AI <^ 1 and in the small-coupling parameter 
C = i/M^ = a^/{/3KM^) < 1, then the CS field equa- 
tions have the solution 



"*Kcrr 



1? 



5 a a cos ( 



\I3M r^ 



1 



12 M 27 AP 

YV ^ 1072" 

2M 1871/2 



sin^ 9dtd4> , 



(28) 



to second order in a/AI and to first order in assum- 
ing no matter sources. These equations employ Boyer- 
Lindquist coordinates {t, r, 9, 0) and ds'^^^.^. is the Kerr 
line element. The solution for may also include an 
arbitrary additive constant, but this constant is unim- 
portant, since only derivatives of d enter the CS field 
equations. Recently, the same solution has been found 
to linear order in a/A4 in the Einstein-Cartan formula- 
tion of the non-dynamical theory [48 1. 



The divergence of the field equations reduce to 
VT^'^' = 0. This is because the divergence of the Ein- 
stein tensor vanishes by the Bianchi identities. Mean- 
while, the divergence of the C-tensor exactly cancels the 
divergence of the CS scalar field stress-energy tensor, 
upon imposition of the equations of motion [Eq. (j25b[) ] . 
Therefore, test-particle motion in dynamical CS gravity 
is exactly geodesicQ. This result automatically implies 
the weak-equivalence principle is satisfied. 

The gravitational perturbation only possesses two in- 
dependent, propagating degrees of freedom or polariza- 
tions. Jackiw and Pi showed that this was the case in 
the non-dynamical theory (32j , while Sopuerta and Yunes 
did the same in the dynamical version [23 . One can also 
show easily that a transverse and approximately traceless 
gauge exists in dynamical CS gravity. The trace of the 
field equations take the interesting form 



(29) 



where R is the Ricci scalar and T is the trace of the stress- 
energy tensor. Notice that the trace of the C-tensor van- 
ishes identically. 

In vacuum (T^'J"' = 0) and when expanding to linear 
order about a Minkowski background, Eq. reduces 
to 



2k 



(30) 



where h = if'''^h^^ is the trace of the metric perturbation, 
□ is the d'Alembertian operator with respect to the back- 
ground metric and § is the background scalar field. Since 
the latter must satisfy the evolution equation [Eq. ()25bp ] , 
we immediately see that cx a/ (3. This means that the 
right-hand side of Eq. (PU]) is proportional to C,. To ze- 
roth order in the small-coupling approximation, h then 
satisfies a free wave equation and can thus be treated as 
vanishing. Deviations from the trace-free condition can 
only arise at C(C) and they are suppressed by factors of 
the curvature tensor, as d must satisfy Eq. (j25b|) . Ap- 
proaching J^+, the right hand side of Eq. (PU)) vanishes. 
This allows one to impose TT gauge at future null infin- 
ity. 



V. EFFECTIVE STRESS-ENERGY IN CS 
GRAVITY 

The perturbed Lagrangian for the Einstein-Hilbert ac- 
tion has already been calculated, so here we need only 
consider the contribution from Sqs- At 0{e^), there are 
a large number of terms generated (we used the package 



* This statement is true only in the absence of spins, since oth- 
erwise the CS effective worldUnc action would contain new self- 
interaction terms. 
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xPert |49l - l53l | to calculate the perturbations). Many of 
these terms are irrelevant when considering their contri- 
bution at J^+. 

Let us classify the types of terms that arise in Scs- -^f 
0{e^), these are of two types: 

1. the second-order part of one field, or 

2. the product of first-order parts of two fields. 

As mentioned in lSec. II Al terms containing the second- 
order part of one field are linear in a short-wavelength 
quantity, which vanishes under averaging. Thus we only 
need to consider the latter case. There are five fields in 
S'cs {V~S' ^1 ^iid *^)' so o^"^ would at first think 

that there are (2) = 10 types of terms arising; however, 
from the definition of the Levi-Civita tensor, we have 

v/^e°^^'' = sign(g)[a/3MH, (31) 
where [aPfj.h'] is the Levi-Civita symbol, which is not a 
spacetime field. The combination \/—g £°'l^t^'^ therefore 
has no perturbation, and there are only three spacetime 
fields which contribute. We are left with only (2) = 3 
possibilities for the types of terms that could appear, cor- 
responding to two perturbed fields and one unperturbed 
one among the set ("i^, RapaS, Rapas)- Two of these pos- 
sibilities are actually the same by exchanging the two 

copies of RaPcrS- 

Therefore, we are left with only the following two types 
of terms in the CS Lagrangian density: 




(32a) 
(32b) 



where W saf) is the first-order perturbation to the Rie- 
mann tensor, in terms of h . Variation of these terms 
with respect to the background metric yields the CS con- 
tributions to the effective stress-energy tensor, while vari- 
ation with respect to the metric perturbation yields CS 
corrections to the first order equations of motion. 



A. Variation with respect to the Perturbation 

Just as in GR, the final expression for the stress-energy 
tensor must be put on-shell by imposing the equations of 
motion. The first-order equations of motion of dynamical 
CS gravity, in vacuum and at are 



Uh 



Al 



(33) 



Imposing these equations of motion is easier when tak- 
ing advantage of the weak-coupling limit, Cgw <C 1, where 
Cgw = § / [nXavi) quantifies the size of the deforma- 
tion away from GR. Let us then expand the metric per- 
turbation in a Taylor series 



n=0 



\SGVi) Hfj.1, ■ 



(34) 



To zeroth-order, it is clear that Eq. ((33|) reduces to 

^kfii, = Oj which is the standard GR equation of mo- 
tion. To next order, the the leading-order piece of the 
right-hand side vanishes and one is then left with 



K 



V Zi,.)(o) 



(35) 



In the remainder of this section, we drop the superscripts 
that indicate Ccw-ordering. 



B. Variation with respect to the Background 

Let us first discuss terms of type >C^^ under variation 
with respect to g^''. From [Sec. II B[ only total derivative 
terms arise from SR'^sap, and these vanish upon aver- 



aging. The remaining terms contain R"' 



SaP 



which must 



vanish at J^+. Thus, as mentioned before, terms in the 
effective action which contain curvature tensors do not 
contribute to the effective stress-energy tensor at 

We are then only left with Cj^j^. Writing these in terms 



of h , the effective action reads 



„tt(2) 



/-off, 1 _|_ /-off, 2 



(36) 



where 



!ff,l _ 



PLla viVpft^^ , (37a) 

^cs' = e"^^^ ^ ^P^pkr ^cr^sh% . (37b) 

Naively, one might think that these expressions lead to 
an effective stress-energy tensor at 0{X^^). This is pre- 
mature, however, as there can be a cancellation of A^^" 
terms that lead to a less steep wavelength dependence. 
One should try to move as many derivatives away from 
the perturbed quantities as possible before procee ding . 
In fact, we know that this must be possible from [32| : 
the Pontryagin density can be written as the divergence 
of a 4-current, so at least one derivative can be moved off 
of h . This automatically implies that there cannot be 
Aqw terms in the effective stress-e nerg y tensor, as shown 
explicitly by Yunes and Sopuerta [g^ . 

Let us transform Ccs^ in the following way. The Levi- 
Civita tensor is contracted onto two derivative operators 
(V^ and V5). One may integrate by parts to move one 
of these derivative operators onto the remaining terms 
in Eq. p7ap . This generates two types of terms: one 
with three derivatives acting on the metric perturbation 
and one with one derivative on the CS scalar (the term 
acting on the Levi-Civita tensor or the determinant of 
the metric vanishes by metric compatibility). Let us fo- 
cus on the former first. Because of the contraction onto 
the Levi-Civita tensor, only the antisymmetric part of 
the second derivative operator would contribute. Such a 
combination is nothing but the commutator of covariant 
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derivatives, which can be written as the Riemann tensor, 
and thus vanishes at J^^ . The remaining term with a 
covariant derivative of the CS scalar does not generically 
vanish. Dropping terms proportional to the Riemann 



tensor, £cs ^ becomes 



/•off, J 
■'-'CS 



(38a) 



Equation (j37bp can be analyzed with the property dis- 
cussed in ISec. II Dl Lorenz gauge may be imposed at 
the level of the action for the purposes of calculating the 
effective stress-energy tensor. This means that if one 
integrates by parts, moving Vo- and Vp onto remain- 
ing terms, the only term that survives is proportional 
to "iJ, as the divergence of Ji^ vanishes (after commut- 
ing derivatives, dropping Riemann terms and imposing 
Lorenz gauge). Thus £cs^ becomes 



VpV^Z? V„A/V^^/ . (38b) 



Let us now evaluate T^!^^^ on shell. Since T^^^^ is 
0((Cgw)"), imposing the equation of motion Eq. (|35|) will 
introduce terms of O^Cgw), which are kept since they are 
the same order as Tcsfiv and Tcsfiiy- We can also impose 
a gauge condition. We have already imposed the Lorenz 
gauge throughout at the level of the action. We may 
further specialize this to the TT gauge. While the TT 
gauge may not be imposed globally, it may be imposed 
at J^^, where the effective stress-energy tensor is being 
evaluated. In TT gauge. 



1_ r ^pr"/^ 



c!!,2 



(42a) 



With these simplified Lagrangian densities at hand, we 
can now compute the total effective stress-energy tensor 
for GWs in CS gravity: 



where 



rjicll 



where T^s^y and T^s^^y are due to the variation of £cs ^ 
and >Ccs'^ respectively. These expressions are 



off, 2 



^ctf.l 



CSflU 1 



(39) 



rj-ioff, 1 



(40a) 



rpc!!,2 _ _ 2 2 
MT/il/ ^ ^ 



P—a^ 



(42b) 
(42c) 



Integrating by parts, imposing the equations of motion 
Eq. (j35p . and integrating by parts again where appro- 
priate, these contributions to the effective stress-energy 
tensor at are 



and 



yoff,2 ^ 2 



C. Imposing the On-Shell Condition 



(40b) 



The equation of motion may be imposed anywhere 
Oh^p may be formed in T^^^^ via integration by parts. 
There is no contraction of derivative operators onto each 



other in Tcspjy, so it remains unchanged. In the final two 
terms of T^s^^ , the derivative operator Vp may be moved 
onto Va")? V hi^x- This would generally make two terms, 
but the term proportional to VqI? '^h^x is O(Cgw) 
ative to the Isaacson piece, so we only keep one term. 
This gives 



rpa{{,l 2 



.Tnoff,2 5a/3.Tncff,l 

^MTflV ^ 6pi/6 ■ 



Finally, we may write an expression for T^^^^ 
after imposing the equations of motion. 



Tncff Tncff 



at 



(42d) 
(42e) 



(43a) 



Tnoff , 1 



■J!,2 



T^ott, 1 



■TcsV., (43b) 



!tf , 1 



.(77 



y)7 



where 5T^!^^^ contains the Chern-Simons correc- 
tion at ©(Cgw)- The summands are taken from 
Eqs. (|40bl) . dtl]), (|42dl) . and P2e)) . Putting them together 
for convenience, the final result is 
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ST"" = 



-6 r a 



(44) 



In the above, we have organized the terms by their 
scaling with powers of wavelength. The first line con- 
tains terms which scale as and A^^; the first of 
these corresponds to a "mass" term in the effective stress- 
energy tensor. Both of these scale more slowly with in- 
verse wavelength than the GR contribution, so they arc 
subdominant. The next three lines have the same scaling 
with inverse wavelength as GR, A^^. The final line scales 
more strongly with inverse wavelength, A^^. This term 
in principle could dominate over the GR term in the high 
frequency limit. 

Notice that the effective stress-energy tensor presented 
here is applicable to both the dynamical and the non- 
dynamical version of GS gravity. Also note that if were 
a constant, rather than a function, the effective stress- 
energy tensor would be identical to that of GR (which 
is expected, since, in that case, the modification to the 
action is purely a boundary or topological term). 



flat, Gartesian spatial slice, appropriate to J^+) must be 
finite. This restricts Vi? to fall off at least faster than 
7.-3/2 then conclude that the GS correction to the 
energy- momentum tensor must vanish at J^'^ , as T^t/'ii/, 
Tmt'^i/j Tcs^i, and T^s^i, decay at least as r"'^/^ or faster. 

The only contribution at ^+ to the effective stress- 
energy of GWs in dynamical GS gravity which decays as 
is the GR part. 

Again, we stress that this only accounts for the out- 
going GW radiation. However, the same argument as 
in ISec. II G] holds; the correction to the energy flux ab- 
sorbed by trapped surfaces is only important at the end 
of an inspiral, both in GR and deformations away from 
GR. This is supported by the small velocity, small mass 
ratio expansion of [l^l (see also [3, |4^ ) . 



D. In dynamical CS gravity 

From asymptotic arguments, we can argue that STl"^^^ 
does not contribute to dissipation laws at in the dy- 
namical version of GS gravity. As mentioned in lSec. II Gl 
the dissipation of energy, linear and angular momentum 
of a system is computed by integrating components of the 
stress-energy tensor on a 2-sphere at J^~^ . Since the area 
of the 2-sphere grows as , for the dissipation integrals 
to be finite, the components of the stress-energy tensor 
must fall-off at least as r~^. In fact, only the part 
of the stress-energy contributes as one takes the r — > 00 
limit. Therefore, any part of the stress-energy tensor that 
decays faster than does not contribute to dissipation 
laws. 

The GS correction to the effective stress-energy tensor, 
'^^cspi/j always falls off faster than in the dynamical 
theory. To see this, we must analyze the behavior of 1?, 
which is restricted. This restriction comes from demand- 
ing that the field sourced by an isolated system and 
in an asymptotically flat space contains a finite amount 
of energy. The energy in 1? is computed by integrating 
the time-time component of T^l^^ on a hypersurface of 
constant time and over all space. For the energy to be 
finite, the integral J°° (V?9)^ r^dr (in an asymptotically 



E. In non-dynamical CS gravity 

In the dynamical theory, since the scalar field must 
carry a finite energy, we were able to argue for the van- 
ishing of (JTcs/ii/ at III fliG non-dynamical theory, 
there is no such demand and no further simplification 
can be made beyond the vanishing of Tllt^ . However, for 
a particular choice of d field, the effective stress-energy 
tensor may be evaluated. We demonstrate this below. 

In the canonical embedding 

The canonical embedding of non-dynamical GS gravity 
is given by (32j 

v^ = V^^^ (1/^,0,0,0) , (46) 

in Gartesian coordinates in the asymptotically-flat part 
of the spacetime. Approaching infinity, this yields 

V„V^i9 = (47) 

so by extension Tcs^f^ = 0, the first-order equation of 
motion becomes D/i^y = + C'(Cq„), the final two terms 
of Tci'f^i/ vanish, and TM-rjuy = = TuT^iy. Notice that 
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here there is no amphtude birefringence in flat spacetime 

as 1? = Hi,!!!,!!!,!!!. 

The first term of Tcs^\, is the only O(Cgw) correction 
which survives. The total stress-energy tensor in the 
canonical embedding of non-dynamical Chern-Simons 
gravity at J'^ ^ with this correction, is 



ST"" =-e^- 

2 « 



2/x 



(48) 



where e*-'*^ is the Levi-Civita tensor on the 3-spacc or- 
thogonal to {d/dt)^ , and the sign change arises from the 
factor of sign(g) in Eq. (PT|) . 



From the form of the correction 5T^l^^^ we can briefly 
mention the leading modification to radiation reaction in 
a binary inspiral at Newtonian order. At this order, there 
is no modification to the trajectories of the two bodies 
from the •& field. Since the first-order equation of motion 
is identical to that of GR at order C(Cgw), the leading 
solution to h^i, is the same as in GR, h^,^ = h^^, . 

Inserting this solution in TT gauge into ST^'^^^ , the en- 
ergy, linear momentum, and angular momentum radiated 
by the system can be computed. Adopting a Cartesian 
coordinate system at asymptotic infinity, the correction 
to the radiated quantities is given by 



(49c) 



where 5Tcs^^l is the part of 6T°"^j^ which decays as 
[5^. In evaluating these integrals, the only angular 
dependence is in factors of rii or Xi. An angular integral 
of an odd number of such factors vanishes, while an inte- 
gral of an even number of them reduces to a symmetrized 
product of Kronecker delta tensors. These factors arise 
explicitly in the definitions of Eqs. and from spatial 
derivatives acting on h^^^ in r°'J. The most important 
difference between T^%^^ and ST°"^j^ is the parity of the 
number of derivatives, which leads to the following be- 
haviour. 

In GR, the leading contribution to E"^^ is from the 
(mass quadrupole)^ combination. Compare this with the 
same integral for ST°"qq, where the (mass quadrupole)^ 
term has an odd number of factors of n^, and thus van- 
ishes. The leading contribution is then from the product 
of the mass quadrupolc and mass octupole. 

The same situation takes place in calculating J^. In 
GR, the leading contribution is from the product of mass 
quadrupole with itself. In the correction from CS gravity, 
the mass quadrupole squared term has an odd number 
of factors of nf, the dominant contribution is again from 
the mass quadrupole times the mass octupole. 



Finally, the situation is different in the calculation of 
Pi. In GR, the quadrupole squared contribution to Pi has 
an odd number of rii factors. The dominant contribution 
is from the mass quadrupole times the mass octupole. 
However, for the CS correction, the quadrupole squared 
term has an even number of factors of n;. Using 



~ TT I 

h = T^.T 



this evaluates to 
5Pf' 



UO'KK'^jJ. 



(50) 



(51) 



where is the reduced quadrupole moment of the mat- 
ter, andzf"^ = {d/dtyiij. 

For a binary in a circular orbit about the S-axis with 
masses toi, TO2, total mass m = mi+7Ti2, symmetric mass 
ratio 77 = TOlTO2/"^^, separation d, and orbital frequency 
w, we find the momentum flux correction to be 



SP^ 



8a 



(52a) 



or, in terms of the velocity v ~ Lod, with Kepler's third 
law = m/d, 



5Pi 



128 



(52b) 



where notice that the quantity in parentheses is di- 
mcnsionless. This is to be compared with the leading 
momentum luminosity in GR, which is proportional to 
P^^ oc rj^v^^Sm/m, where 5m = mi — 7712 [53]. Although 
the GR effect is two powers of v stronger, it depends on 
the difference in masses, whereas the non-dynamical CS 
correction only depends on the total mass. This implies 
that in the limit of comparable masses mi m2, the re- 
coil velocity would not asymptote to zero in CS gravity, 
as it does in GR for non-spinning binaries. 

A physical interpretation of this effect is related to 
the parity-violating nature of the theory. When one 
chooses a canonical embedding, the action becomes 
parity-violating as the Pontryagin density is parity odd. 
The embedding coordinate chooses a (temporal) direc- 
tion to which the modification to the Einstein equations 
can couple to, inducing a new term in the stress-energy 
that is proportional to the curl of the metric perturba- 
tion. Because kicks are predominantly generated during 
merger, the CS modification is indeed dominant over the 
GR result, leading to the first, non-linear, strong-field 
modification computed in CS gravity. 



VI. EFFECTIVE STRESS-ENERGY TENSOR 
OF MODIFIED GRAVITY THEORIES 

Let us now consider a broader class of modified gravity 
theories. There is an infinite variety of GR modifications 
one could construct. However, there are several proper- 
ties that are desirable and that we require here: 
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1. Metric theories: the action depends on a symmetric 
metric tensor that controls the spacetime dynamics. 

2. Deformations of GR: analytically controllable and 
small corrections to the Einstein-Hilbert action 
with a continuous GR limit. 

3. High-Rank Curvature: corrections depend on 
quadratic or higher products of the Riemann ten- 
sor, Ricci tensor, or Ricci scalar. 

4. Minkowski stable: the theory must admit 
Minkowski spacetime as a stable vacuum solution, 
and future null infinity should be asymptotically 
flat for isolated matter spacetimes. 

Besides the metric, there may be new fields introduced 
which are considered part of the "gravity sector" . This 
distinction means that said fields are not minimally cou- 
pled, i.e. they may be coupled to connection and cur- 
vature quantities. These additional fields may be of any 
spin: scalars, spinors, vectors, etc. For simplicity, we will 
only consider scalar fields here, but the results may also 
be extended to higher spin fields. Scalar fields are well- 
motivated from quantum completions of GR, e.g. moduli 
fields are common appearances in string theoretical mod- 
els M. 



A. Action 

In defining a modified gravity theory, let us consider 
what terms may arise in the action. These terms must in- 
clude the Einstein-Hilbert and matter terms, along with 
modifications built from additional scalar fields and cur- 
vature invariants. Additionally, it ought to contain a dy- 
namical term for the scalars that couple to the curvature 
invariants, as we will motivate in Sec. IVIBI 

In principle, there are an infinite num- 
ber of curvature invariants to consider. The 
first few of these are simple to construct: 
A, i?, i?^, Vpi?V^i?, R^uR^'^ , RapiiuR"'^^'^ 1 ■ ■ 
where A is any scalar constant, e.g. the cosmological 
constant. These may be specified by their rank, r, 
which is the number of curvature tensors which are 
contracted together, and by further specifying a list of r 
non-negative integers {Ai, . . . , A^}, where Ai specifies the 
number of derivatives acting on the i"^ curvature tensor. 
For a rank r and case {\iYi=ii there are a finite number 
of independent curvature invariants corresponding to the 
number of ways to contract indices. Thus all curvature 
invariants may be countably enumerated, assigning some 
number n to each independent invariant. 

We here consider only combinations of algebraic cur- 
vature invariants, i.e. Ai = for all cases. This means we 
do not allow modifications that depend on derivatives of 
curvature tensors. Such a simplification is a good one, 
from the standpoint that it automatically guarantees the 
field equations to be no higher than fourth-order. 



Consider then a modified gravity theory defined by the 
action 

S = Seh + 'S'mat + >5'i„t + , (53a) 
where Si) is the canonical kinetic term for 



13 



d^A 



(53b) 

with V an arbitrary potential function; and where S'int 
is the interaction term between the scalar and some 
algebraic combination of curvature tensors, for example 



"5'int , 1 

5'i„t,2 



a2 / d^x 




or generally 



(53c) 
(53d) 
(53e) 

(53f) 



with / an arbitrary "coupling function" and TZ an alge- 
braic combination of curvature invariants. Alternatively, 
notice that we could have assigned each term propor- 
tional to ai a separate coupling with its associated 
kinetic and potential terms. The arguments presented 
below would also hold for such constructions. 



B. Dynamical scalar fields 

The requirement for the scalar d to be dynamical arises 
from demanding diffeomorphism invariance in the theory. 
Consider the infinitesimal transformation of the action 
under a diffeomorphism generated by the vector field v'^ . 
Specifically, look at the terms containing -d, i.e. the sum 
•Smod — Sint + S§ . The infinitesimal transformation under 
the diffeomorphism is 



6^ 



where is the interaction Lagrangian density, is 
the kinetic Lagrangian density, and Cv stands for the Lie 
derivative along w^. 

For a theory to be diffeomorphism invariant, the in- 
finitesimal transformation in the total action must van- 
ish, dS = 0. Since £1,1!) may be arbitrary for some ■& and 
some v^^, the functional multiplying C^'d must vanish for 
dS to vanish. This means 



6§ 



(£.„, + £^) = . 



(55) 
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When the scalar field has dynamics, i.e. /3 7^ 0, then 
Eq. ([55)1 is identical to the equations of motion of the 
field 1? and is therefore automatically satisfied. However, 
if the field is not dynamical, /3 = 0, then Eq. ([55)1 gives 



/'(^)7^ = o. 



(56) 



Except in the case where /'{"d) — 0, this is an additional 
constraint on the geometry of spacetime, namely that 
TZ = 0. Given that the equations of motion already satu- 
rate the number of equations for the degrees of freedom 
present, this would be an ovcrconstraincd system. This 
is in fact the case in the non-dynamical version of CS 
gravity, as discussed in [13, Ull . We therefore only ad- 
mit dynamical scalar fields, or terms with no scalar field 
dependence (/(i^) = const.). 



C. Special cases: zeroth and first rank 

Before doing a calculation for a general curvature in- 
variant TZ, let us briefly discuss some special cases. As 
we will see, curvature invariants of zeroth and first rank 
will not be considered. 



Zeroth rank 

At zeroth rank, there is only one algebraic curvature 
invariant: a constant. The non-constant part of /o(i^) 
may simply be reabsorbed into the potential Vl^). This 
gives a minimally coupled scalar field, which may be ab- 
sorbed into S'niaf The constant part of /o leads to a 
"cosmological constant." Since we are only considering 
theories which are Minkowski stable and asymptotically 
flat, this cosmological part must vanish. 



First rank 

There is only one algebraic curvature invariant of rank 
1, the Ricci scalar R. If we allow / to be a non-constant 
function, we would have a classical scalar-tensor the- 
ory, akin to Brans-Dicke theory (see e.g. [l|). The ef- 
fective GW stress-energy tensor for scalar-tensor theory 
has been computed, for example in Brans-Dicke theory 
(see [i^l), so we do not consider it here. Since we already 
include the Einstein-Hilbert term in Eq. (|53ap , there can 
be no additional term linear in R without affecting New- 
ton's constant. Thus we only consider quadratic and 
higher rank curvature invariants. 



D. Cubic and Higher Ranks 

At cubic rank, one can easily show that there 
are five algebraic invariants that may not be 
factored as products of lower rank invariants 



{R^ifR pR^^: R^ R ^ Rpai^p 7 R ^ ^vR^ paR^ q/3, 
*RP''^"'R'^^^yRpa^\, and * Rp^'^f^ RP o^R" p) and four that 
may be factorized (i?^, RRp,,R^"', RR^^^pR^""^ , and 
R*RP^P''R'' p^i,). The arguments that follow work for 
all of them, so for concreteness we choose just one: 
RP" ^R'^ pRP p,. The modification to the action arising 
from this term is 



^ f{^)R^',R%RPp 



(57) 



The contribution from this term to the effective action 
at second order is 



^off(2) 



f'{d)d R^'^R^RP^ 



-f 3/(^9) R^,R%RPp]+if'{d)d R^^.R^pRPp 



+3f{d)R'',R%RPp\ . (58) 

Immediately we see that all terms have at least one power 
of background curvature tensors. This means that each 
term can be written similarly to an earlier example in 
Sec. Ill C[ in Eq. pT|) . When evaluating this effective 
stress- energy tensor at , all of the background curva- 
ture tensors vanish. This automatically implies that cu- 
bic and higher rank terms in the action do not contribute 
to the effective stress energy tensor at asymptotically- 
flat, future null infinity. 

The stress-energy tensor is then given by the 
MacCallum-Taub tensor, Eq. (j20bp . which need not be 
identical to the GR one yet, as one must first impose the 
first-order equations of motion at J^"*". These equations 
could be modified by the introduction of higher-order op- 
erators in the action. Let us analyze such equations again 
through the example of Eq. ([57)) . As the calculation de- 
pends only on the rank of the curvature invariant ap- 
pearing in the action and not on its specific form, the 
results shown below extend to all cubic and higher rank 
algebraic curvature invariants as well. 

The equation of motion arising from Eq. ()57|) is 



KGp, + Zaf{d)RppRf'^R\ - -g^Ji^)R''pR^^R\ 

-2V^V(^ (/(^?)i^\)i?^)] = r;r + t^V • (59) 

The important feature to note is that all terms containing 
a, that is, all terms deforming away from GR, arc cubic or 
quadratic in curvature tensors. This is a general feature: 
from a term in the action of rank r, terms in the equations 
of motion will be of rank r and rank r — 1. 

Now consider evaluating the first-order equations of 
motion at asymptotically-flat, future null infinity, which 
we need in order to put the MacCallum-Taub stress- 
energy tensor on-shell. We will not write out the full 
first-order equations of motion; it suffices to say that the 
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modification terms (those terms containing a) are of rank 
r, r — 1, and r — 2 in tlie first-order equations of motion. 
When going to only the terms of rank survive, 

e-g- Oh^^. 

Immediately we see that the only modifications to the 
action that affect the first-order equations of motion at 
are those of rank 2 and lower. Thus for modifications 
that are cubic and higher, the first-order equations of 
motion at are simply those of GR, QA^jy = 0. 

Inserting this asymptotically-flat, on-shell condition 
into the MacCallum-Taub stress-energy tensor yields the 
Isaacson stress-energy tensor. Cubic and higher rank 
modifications to the Lagrangian do not modify the ef- 
fective stress-energy tensor due to GWs. We again em- 
phasize that radiation reaction will still be different in 
a higher order theory because of different motion in the 
strong field, additional energy carried in the scalar field 
1?, and energy carried down horizons being different. But 
the energy of a GW at ^+ is the same as in GR. 

E. Quadratic terms 

Let us now consider quadratic deformations to the ac- 
tion, as these are the only ones left to study, and let us 
classify the types of modifications possible. There are two 
important characteristics that we use for such a classifi- 
cation. The first depends on the nature of the curvature 
quantity TZ. This quantity may either be topological or 
not. A curvature quantity that is topological may be ex- 
pressed as the divergence of a current, TZ = V^/C''. As 
we mentioned in Sec. IIV[ Eq. p4)). the Pontryagin den- 
sity, is a topological curvature invariant. In metric 
gravity, the only other non- vanishing, algebraic, second 
rank curvature invariant that is topological is the Gauss- 
Bonnet term, 

g = Ro.Pf.uR''^^'' - 4i?^.i?'"' + i?" , (60) 

as the Nieh-Yan invariant vanishes in torsion-free theo- 
ries. 

The second characteristic we can use to classify theo- 
ries is the behaviour of the scalar field i?, which depends 
on the potential V{d). The two possibilities are a po- 
tential that is flat, V{d) = 0, or one that is non-flat, 
V varying with A flat potential does not choose out 
any preferred values of the scalar field, whereas a non- 
flat potential must be bounded from below for stability, 
and thus has a global minimum (or several minima) . The 
presence or absence of a preferred field value is important 
in the limit going to 

For non-flat potentials, without loss of generality, the 
global minimum can be shifted to = by simulta- 
neously shifting the potential function and the coupling 
function /. Such a shift does not affect the derivative 
term in the kinetic term for iJ, since it simply adds a 
global constant to the field. The only problematic situ- 
ation is if the global minimum is in the limit ?9 — >■ ±oo, 
which we do not allow here. 



We begin by discussing the asymptotic behaviour of "d, 
which satisfies the sourced wave equation 

li{Ud~V'{d))^-af'{d)'R. (61a) 

At the right hand side vanishes. Furthermore, if 

we are interested in static or quasistatic background so- 
lutions for around which we can expand, time deriva- 
tives in the d'Alembertian will vanish leaving only the 
Laplacian, 

V^i?-y'(i?) =0. (61b) 

For a non-flat potential, V^d) ^ 0, the zcroth order 
asymptotic solution will be i9 going to the minimum of 
the potential, which we have shifted to = 0. 

For a flat potential, the background equation of motion 
for d, Eq. (|61bp . at ^+ becomes 

V^i? = . (61c) 

There are two asymptotic solutions: -d asymptotes to a 
constant or d asymptotes to a function linear in Carte- 
sian coordinates. The latter case would contribute a con- 
stant stress-energy tensor at . This would lead 
to an asymptotically de Sitter spacetime, not an asymp- 
totically flat spacetime. Therefore, we only consider the 
case where asymptotes to a constant. 

The equation of motion Eq. (j61cp does not determine 
to what value 'd asymptotes. A boundary condition is re- 
quired in this case. Again, without loss of generality, for 
some given asymptotic value determined by some bound- 
ary condition, the field and coupling function f{'&) may 
be shifted so as to redefine the asymptotic value to be 
I? — >■ without changing the physics. 

A boundary condition is not required if the theory is 
"shift symmetric." In a shift symmetric theory, the trans- 
lation operation i? — > -f- c, where c is a constant, leaves 
the equations of motion invariant. Such a theory, there- 
fore, must have equations of motion that depend only on 
the derivative V^i?. Such is the case, for example, if the 
action depends on a topological term multiplied by some 
scalar field, /(i9)V^/C^, as then the action can be rewrit- 
ten as (V^/(7?))/C^ via integration by parts. Of course, 
in this case, the potential must also be flat and / must be 
linear for the theory to be shift-symmetric. Such types 
of corrections arise naturally in the low-energy limit of 
string theory ^UMl. 

Let us rewrite the action and split the interaction term 
into a dynamical and non-dynamical part. Since we can 
always shift the field, potential, and coupling function so 
that the asymptotic value is i? — 0, let us define 

a' = af{0) (62a) 
F(i?) = fi^) ~ /(O) . (62b) 

Then, the interaction term in Eq. (|53f1) may be rewritten 

as 

= 5„-d + 5,,„ (62c) 
= a' I d^a;^/^ TZ + a [ d^xy^ F{i})TZ . 
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Field 


Asymptotic form 


Field 


Asymptotic form 






V (/ 




r 






At least r"^/^ 


R 






At least r"^/^ 
















-1 

r 











TABLE I. The asymptotic forms of fields appearing in the 
effective action for a rank 2 modification to the action. All 
tensor indices have been suppressed. 



The first term is the non-dynamical part, i.e. the part 
that does not couple to the scalar field, while the second 
part is the dynamical part. If 7?, is a topological curva- 
ture invariant, then the first term in Eq. (|62cp does not 
contribute to the equations of motion, as it is the integral 
of a total derivative. 



1. Dynamical contribution 

Let us perturb S^yn to second order to calculate the 
contribution to the effective action, keeping in mind that 

d and h do not contribute. This part of the effective 
action is 



?«ff(2) 

^dyn 



^-F{^) + F'{d) d^n 



+ ^ {h F'{&)^ + F"{d) Te, 



(63) 



To determine the contribution to the effective stress en- 
ergy tensor at ^+ , again analyze the asymptotic form of 
all of the fields appearing in Eq. The asymptotic 

forms are summarized in Table ID 

The simplest way to see that the dynamical part of 
the effective action docs not contribute at J^"*" is to ex- 
amine the asymptotics of TZ, TZ, and TZ. Since curvature 
tensors Ra/j^u SuTC tidal tensors, they goes as r~^; since 
TZ contains two curvature tensors, it scales as . The 
slowest decaying (i.e. leading) part of TZ roughly comes 
from RR (with indices suppressed); the leading part of R 
is V h, which, being radiative, goes as r~'^ . This means 
that IZ ^ r^*. Similarly, the leading part of TZ goes as 
(V^/i)(V^/i), so7^-r-2. 

Examining the effective action, Eq. (|63p . we see that 
there arc no terms that decay as which are the 

only ones that can contribute to the GW effective stress- 
energy tensor at .y^. Any term with TZ ot: TZ already 

decay too quickly; only terms with TZ could remain, and 



only if they were multiplied by terms that asymptote as 

o{i). ^ 

Having performed the splitting into the dynamical and 

non-dynamical parts, TZ is multiplied by ^~gF{'d) in 
the effective action. This splitting was specifically con- 
structed so that F(0) = 0. Since F must be differentiable 
at I? = 0, J^(i?) must go to zero at least as fast as i) goes 
to zero, which is at least r~'^/'^. 

We have thus shown that the dynamical part of the 
interaction term does not contributed to the effective 
stress-energy tensor at directy. However, it could 
still contribute indirectly through the imposition of the 
first order field equations. We examine this in a later 
section. 



2. Non-dynamical contribution 



Let us now consider S'„_d in Eq. (|62cp . This term gener- 
ically contributes to the effective stress-energy tensor of 
GWs at To show this contribution, consider the 

general rank 2 modification as the linear combination of 
the four independent rank 2 curvature invariants 

aTZ = aiR^+ aai^^^.i?^'' + a^Rapt^^R"''^^'' + a4*RR, 

(64a) 

and absorb /(O) into the coefficients a'^ in the non- 
dynamical part, 

a'TZ EE a[R^ + a'^Rf.^R^"' + a^Ra(i^,uR°'^^'' + a4*RR. 

(64b) 

Note that this form also includes the Weyl squared invari- 
ant, which is a dependent linear combination of the above 
terms, C"^''^C„;3^, = R^i - 2R'^PR^p + R'^f^'^" Rc^fi^., 
which is considered in (26l . [27j . The Pontryagin den- 
sity *RR, being a topological invariant, does not con- 
tribute to the action in 5'„_d, so we may drop the final 
term. Similarly, if the linear combination is propor- 
tional to the Gauss-Bonnet (or Euler) invariant, which 
has ai = 1 = a^, a2 = —4, then TZ would be topological 
and there would be no contribution to 5'„_d and hence no 
contribution to the effective stress-energy tensor of GWs. 

The calculation of the effective action for the non- 
dynamical term is straightforward but long, so we do not 
show the steps here. An outline of the calculation is to 
perturb ^J—gTZ to second order; the only parts that may 
contribute to an effective stress-energy tensor at ^+ are 
of the form y/—gRR, where again we have suppressed in- 
dices on the perturbed curvature tensor R. This is calcu- 
lated in terms of the trace-reversed metric perturbation 
Ji^ . As before, Lorenz gauge may be imposed at the level 
of the action. All terms that remain will be of the form 
VaV^^^'' VkVaA'"^ with all indices contracted to form 

~ flu 

a scalar. If any derivative is contracted onto h , by in- 
tegrating by parts and commuting covariant derivatives, 
one may form the Lorenz gauge condition V^/i — and 
ignore the term in the effective action. Thus the only sur- 
viving terms have derivatives contracted together, which 
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— ~ /it/ — ~ _ ~ _ ~ 

can be put into one of two forms, Oh Oh^^^ and Oh Oh. 

After the cxplieit calculation, the prefactors are found 

and 



and the trace of this equation is 

nOh ^ 2 (Sai + a2 + as) f{d)OOli . 



(68b) 



ff(2) 



a'^)Oh Oh 



a. 



Oh 



(65) 



Note again that does not appear, and if ai, 02, as 
are in the Gauss-Bonnet ratio, then the effective action 
of Eq. (p5|) vanishes. 

Putting all indices in their natural positions, so as to 
expose implicit metric dependence, and varying the ef- 
fective action of Eq. (|55|) with respect to g'"^ , the contri- 
bution to the effective stress-energy tensor is 



(a'l - a'3) Oh (Oh^^ - V^V.h 



(66) 



We then find that the effective stress-energy tensor 
of GWs at J^^ is given by the MacCallum-Taub stress- 
energy tensor (coming from the Einstein-Hilbert action) 
plus the direct contribution from the non-dynamical part 
of the rank 2 interaction term, r°"pj,, 



rpcii rricff 



(67) 



The only remaining part of the calculation is to put the 
stress-energy tensor on-shell, that is, to impose the first- 
order equations of motion at 



3. First-order equations of motion 

We need the first-order equations of motion at J^^ of 
the full theory, including both the dynamical and non- 
dynamical terms of the action. At J^^, however, •& and 
V^i? decay at least as r~^/^ and r^^/^. The only remain- 
ing dependence on d is through f{d) /(O). 

The general zeroth and first-order equations of motion 
are quite long, so we do not reproduce them here, but 
they do simplify as r -> 00. These equations are lin- 
ear in h^^ and S, which are both radiative and decay 
as This r^^ scaling is the asymptotic scaling of 

the first-order equations of motion, as there are terms 
of the form Oh^^^ that appear with no curvature tensors 
or background scalar field 1} multiplying them. However, 
all terms containing have curvature tensors multiplying 
them, so they decay faster than the leading behaviour of 



Again we see that if the a coefficients are in the Gauss- 
Bonnet ratio, the GR equation of motion is recovered at 

This wave equation can be seen to be a massive wave 
equation for the auxiliary variable f^,y = D/i^^, with 
mass m ^ 1/A, where ^ |Q;i|/(0)/K. In the weak 
coupling limit, A/Agw ^ 1, the equations simplify con- 
siderably. This simplification comes from treating the 
solution to the full theory as a deformation away from 
GR; this means expanding the fields as power series in a 
small parameter, namely ( = (A/Aqw)^- As in Eq. (|34p . 
we impose 



00 
71=0 



(69) 



and similarly for other fields, where the zeroth field h^^^ 
is the GR solution. Inserting this expansion in the first- 
order equation of motion Eq. (j68p and matching order by 
order gives 

KOh%'''^ = - {2a, + ^2 + 2^3) 

x(V^V,Q-g^,DQ)A^"^ 

- {a2 + 4a3) f{^)aOh}^J , (70a) 



for all orders n > 0, and 



0, 



(70b) 



for the GR solution. Substituting Eq. (j70bp into 
Eq. (|70ap and iteratively solving the field equations one 
order at a time, we find at all orders that 



0. 



(70c) 



This is the GR first-order equation of motion, and just 
as in GR, we may specialize the Lorenz gauge to the 
TT gauge at This expansion has discontinuously 

turned the massive wave equation into a massless one by 
killing the massive modes in the limit of m — > 00. Such 
an order-reduction procedure, where certain solutions are 
eliminated through perturbative constraints, has been 
shown to select the physically correct ones in all stud- 
ied cases [31I, . 

We can now evaluate the complete effective stress- 
energy tensor of GWs at J^~^. As shown in Sec. IVIE 11 



Keeping only the terms that go as r^^, the first-order 
equation of motion in Lorenz gauge at asymptotically- 
flat, future null infinity is 



iOh 



(2ai + a2 + 2^3) fW (V,,V.D - g^.DD) h 
(a2+4a3)/(^)QQ^^,, (68a) 



To prove that the TT gauge exists at for this theory, the 
proof in Appendix A of Flanagan and Hughes (FH) l64i must 
be extended. Their Eq. (A. 12) must be replaced by our Il68al l 
and a small coupling expansion performed. The result will again 
be I ITOcll . which is identical to FH's Eq. (A. 12). 
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there is no direct contribution from the dynamical part 
of the interaction term. Section IVIE2I showed that the 
non-dynamical part does contribute directly, but impos- 
ing Eq. (|70c[) forces this contribution to also vanish. Since 
the MacCallum-Taub tensor on-shell is equal to the Isaac- 
son tensor, we then have 



rpcff irjicU 

flu \ MTflL' 



n-d/J-Z^ J 



T°" (71) 



That is, the effective GW stress-energy tensor is identical 
to the Isaacson one at for this wide class of modified 
gravity theories. 



VII. CONCLUSIONS 

We have here addressed the energy content of GWs 
in a wide class of modified gravity theories. We focused 
on theories that are weak deformations away from GR 
and calculated the effective stress-energy tensor where 
GWs are extracted: in the asymptotically-flat region of 
spacctime. 

The main calculation tool we employed was the per- 
turbed Lagrangian approach. We demonstrated the cal- 
culation explicitly for GR, recovering the Isaacson effec- 
tive stress-energy tensor. We also explicitly calculated 
this effective tensor in dynamical modified CS gravity, 
where again the result at J^"'' reduces to the Isaacson 
tensor. The features of CS gravity that lead to the ef- 
fective stress-energy tensor being identical to the one in 
GR are the dynamical nature of the scalar field and the 
topological nature of the curvature correction to the ac- 
tion. 

We then generalized this finding to all action modifica- 
tions of a similar nature: a dynamical scalar field coupled 
to a scalar curvature invariant of rank 2 or higher in a 
spacetime that is asymptotically flat. For scalar curva- 
ture invariants of rank 3 or higher, we showed that there 
is no modification to the stress-energy tensor or the equa- 
tions of motion at For rank 2, we calculated the 
contribution to the effective stress-energy tensor and to 
the first-order equation of motion. In the weak coupling 
limit, the only solutions to the first-order equations of 
motion satisfy the GR first-order equations of motion at 
namely Oh^^, ~ 0. Evaluating the effective stress- 
energy tensor on-shell with these solutions leads, again, 
to the Isaacson stress-energy tensor. 



A few caveats are in order. As we have stressed be- 
fore, this result is evaluated at asymptotically-flat, fu- 
ture null infinity, so it does not apply to cosmological 
spacetimes, e.g. de Sitter spacetime. Not all of the en- 
ergy that is lost by a system is carried away by GWs to 
^+: there is also radiation in the scalar field (which is 

calculated straightforwardly from Tj^V), and both GWs 
and the scalar field radiation are lost to trapped surfaces. 
All of these effects must be accounted for in calculating 
the radiation-reaction of a system. Finally, we did not 
address modifications to the action of the form f{'d)R, 
which reduce to a classical scalar-tensor theory. 

There are several avenues open for future work. Con- 
sidering classical scalar-tensor modifications is one pos- 
sible extension. The work should also be extended to 
the next simplest spacetimes, those that are asymptoti- 
cally de Sitter. This is appropriate for calculating GWs 
from inflation, for example. Extending this approach to 
calculating energy lost to trapped surfaces is another pos- 
sibility. 

The most natural application of this work is in tests 
of GR with pulsar binaries and with GWs emitted by 
EMRIs. The former problem requires performing a post- 
Keplerian expansion of the motion of bodies orbiting each 
other. The latter requires knowing the BH spacetime 
(background) solution in the class of modified gravity 
theories and the geodesic or non-geodesic motion on that 
spacctime. Both of these programs require knowledge of 
radiation-reaction in GWs at J^~^, which we have here 
computed for a large class of modified gravity theories. 
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